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Vector operators based on robust order statistics have proved successful 
in digital multichannel imaging applications, particularly color image filter- 
ing and enhancement, in dealing with impulsive noise while preserving edges 
and fine image details. These operators often have very high computational 
requirements which limits their use in time-critical applications. This paper 
introduces techniques to speed up vector filters using the minimax approxima- 
tion theory. Extensive experiments on a large and diverse set of color images 
show that proposed approximations achieve an excellent balance among ease 
of implementation, accuracy, and computational speed. 
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1. Introduction 

Image noise filtering - the process of estimating the original image information from noisy 
data - is a common preprocessing step in image processing and analysis applications, as 
the presence of noise in images not only lowers their perceptual quality, but also makes 
subsequent tasks such as edge detection and segmentation more difficult pp. With the recent 
shift from traditional grayscale imaging to color imaging, numerous filters have been proposed 
for removing noise from color images. An extensive overview of color image filtering solutions 
and their applications can be found in [2] , with detailed performance analysis presented in [3] . 
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An important class of filters for noise reduction in color images is the one based on robust 
vector order statistics [HE]. A typical natural image exhibits strong correlation between its 
red, green, and blue color channels; therefore, treating the pixels of the image as vectors 
avoids color shifts and artifacts in the output of the filter. Since images are nonstationary 
due to the presence of edges as well as noise and blur introduced during the image formation, 
vector filters usually operate on pixels inside a supporting window that slides over the image. 
Desired noise filtering characteristics can be obtained by using vectors with certain ranks 
in the ordered set of pixel values inside the supporting window, as an ordering operation 
performed according to a distance or similarity criterion distinguishes outliers from noise- 
free samples [I]. 

Many researchers have noted the high computational requirements of order-statistics based 
vector filters; however, relatively few studies j6j[T] have focused on alleviating this problem. 
Furthermore, the scope of these studies is limited to the vector median filter [8], which 
has been considered as the gold standard of performance in color image filtering due to its 
robustness and excellent impulsive noise suppression capability [9]. 

This paper introduces techniques to speed up popular vector filters which use vector 
ordering criteria other than the Euclidean distance. In particular, the filtering solutions 
from [101412] involve, respectively, computationally expensive inverse cosine, exponential, 
and logarithmic functions that are evaluated during the filtering process typically millions of 
times. In order to allow the use such filters in time-critical imaging applications, we utilize 
the minimax approximation theory to substitute the abovementioned elementary functions 
with computationally efficient polynomials. Extensive experiments on a large and diverse 
image set show that the presented approximations achieve an excellent balance among ease 
of implementation, accuracy, and computational speed. 

The rest of the paper is organized as follows. Section [2] gives background on minimax 
approximation theory. Section [3] introduces the use of the minimax approximation theory 
in speeding up order-statistics based vector filters. Motivation and design characteristics are 
discussed in detail. Section H] describes the image set, noise models, filtering performance 
criteria, and the experimental setup. Finally, conclusions are given in Section [5j 

2. Overview of Minimax Approximation Theory 

Given a function /, we would like to approximate it by another function g so that the error (e) 
between them over a given interval is arbitrarily small. The existence of such approximations 
is stated by the following theorem: 

Theorem 2.1 (Weierstrass) Let f be a continuous real-valued function defined on [a,b], 
i.e. f G C[a,b]. Then We > there exists a polynomial P such that \\f — P\\ < e, i.e. 
Vze [a, 6], \f(z)-P(z)\<e. 
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This is commonly known as the minimax approximation to a function. It differs from 
other methods, e.g. least squares approximations, in that it minimizes the maximum error 
(e) rather than the average error: 

e=max\f(z)-P(z)\. (1) 

ze[a,b] 

A similar theorem establishes the existence of a rational variant of this method [13]. Let 
n > be a natural number and let 

P n ([a,b]) = {av + aiz+ ... +a n z n : z G [a,b], ajGl, z = 0,l,...,n} (2) 

be the set of all polynomials of degree less than or equal to n. The set of irreducible rational 
functions, R^ ([a, b}), is defined as 

R n m ([a,b}) = !^: p(z) eP n ([a,b}), q(z) G P m ([a, &])} (3) 

where p and q have no common factors. Then [T3] : 

Theorem 2.2 For each function f G C[a,b], there exists at least one best rational approxi- 
mation from the class R 1 ^ ([a, &]). 

This theorem states the existence of a rational approximation r* G RJ^ ([a, b]) to a function 
/ G C[a, 6] that is optimal in the Chebyshev sense: 

max |/(z)-r*(z)|= dirt (4) 

z6[a,6] 

where dist (/, -R") denotes the distance between / and _R^([a, 6]) with respect to some 
norm, in our case, the Chebyshev (maximum) norm. Regarding the choice between a poly- 
nomial and a rational approximant, it can be said that certain functions can be approxi- 
mated more accurately by rationals than by polynomials. Jean-Michel Muller explains this 
phenomenon as follows "It seems quite difficult to predict if a given function will be much 
better approximated by rational functions than by polynomials. It makes sense to think that 
functions that have a behavior that is 'highly nonpolynomial' (finite limits at ±oo, poles, 
infinite derivatives, . . .) will be poorly approximated by polynomials." [14] . 

In this study the Remez Exchange Algorithm, an iterative method that uses Lagrangian 
interpolation to systematically minimize the maximum absolute difference between the given 
function and its polynomial approximation, was used to calculate the polynomials. The reader 
is referred to [T3JIT3] for more information on the minimax approximation theory and [15] 
for the implementation details of the Remez algorithm. 
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3. Proposed Implementations of Vector Filters 

Consider anMxiV red-green-blue (RGB) input image X that represents a two-dimensional 
array of three-component vectors x(r, c) = [xi(r,c),X2(r,c),x^(r,c)] occupying the spatial 
location (r, c), with the row and column indices r = {1, . . . , M} and c = {1, . . . , N}, re- 
spectively. In the pixel x(r, c), the Xk(r,c) values denote the red {k = 1), green {k = 2), 
and blue (k = 3) components. In order to isolate small image regions, each of which can be 
treated as stationary, and reduce processing errors by operating in such a localized area of 
the input image, an y/n x y/n supporting window W(r, c) centered on pixel x(r, c) is used. 
The window slides over the entire image X in a raster fashion and the procedure replaces 
the input vector x(r, c) with the output vector y(r, c) = F(W(r, c)) of a filter function F(-) 
that operates over the samples inside W(r,c). Repeating the procedure for each pair (r, c), 
with r = {1, . . . , M} and c = {1, . . . , N}, produces the output vector y(r, c) of the M x N 
filtered image Y. For notational simplicity, the input vectors inside W(r, c) are re-indexed 
as a set, i.e. W(r,c) = {x^ : i = 1,. . . ,n}(see Figure [Q, as commonly seen in the related 
literature [2}[3]. In this notation, the center pixel in W is given by xc = X( n+1 )/ 2 and in the 
vector Xj = [xn, x i2 , x i3 ] with components x^, the i G {l,...,n} and k G {1,2,3} indices 
denote the position of the vector inside the window and the color channel, respectively. 
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Fig. 1. Indexing convention inside a 3 x 3 window 
3. A. Vector Directional Filters 

The vector directional filter (VDF) family pi)] operates on the direction of the input vectors 
with the aim of eliminating the vectors with atypical directions. This family utilizes the angle 
between the input vectors to order the vectors inside the supporting window. For example, 
the output of the basic vector directional filter (BVDF), the most well-known member of the 
VDF class, is the input vector inside the supporting window whose direction is the maximum 
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likelihood estimate of the directions of the input vectors 



y(r,c) = argmin A^Xj) 

^eW(r,c) U=i I (5) 



A (Xi, x,) = arccos ( 

J \ Il x i|l2ll x j'll2 / 



where A(x.i,Xj) denotes the angle between the two input vectors Xj and Xj and || . || 2 is the 
L2 (Euclidean) norm. Note that in addition to BVDF, the angular function A(., .) was used 
in the design of a number of other filters including the generalized VDF [17], directional 
distance filter [18], hybrid vector filters [19], weighted VDFs [20], data-adaptive VDFs [21], 
and switching VDFs |22j . 

The computational requirements of these filters can be reduced by speeding up the inverse 
cosine (ARCCOS) function, whose argument falls into the interval [0,1] (see Figure |2]). 
Unfortunately, approximating the ARCCOS function in this interval is not easy because of 
its behavior near 1. This can be circumvented using the following numerically more stable 
identity for z > 0.5: 

arccos (z) = 2 arcsin ( a/0.5(1 — z) ) (6) 



where the inverse sine function (ARCSIN) receives its arguments from the interval [0, 0.5] 
(see Figure [3]). Instead of plugging the value of a/0.5(1 — z) into a minimax approximation 
for the ARCSIN function and then multiplying the result by 2, two multiplication operations 
can be avoided if the following function is approximated: 



T = (7) 
arccos (z) = 2 arcsin (r/y/2) 



where the argument r falls into the interval [0, ij a/2] . 

Table [1] lists the coefficients of the fourth degree minimax polynomials that approximate 
the ARCSIN and ARCCOS functions. Since both functions exhibit strong linearity in their 
respective intervals, they can be accurately approximated by polynomials, as indicated by 
the small error values listed in the table. 



Table 1. Fourth degree minimax polynomials for the ARCSIN and ARCCOS 
functions 



Function 


£ 


ao 


01 


12 


03 


C14 


ARCSIN 


2.097814e-05 


2.097797e-05 


1.412840 


1.429881c-02 


6.704361e-02 


6.909677c-02 


ARCCOS 


1.048949e-05 


1.570786 


-9.990285c-01 


-1.429899e-02 


-9.481335e-02 


-1.381942c-01 
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Fig. 2. Function arccos(z) in the interval [0, 1] 




Fig. 3. Function arcsin(z) in the interval [0, 0.5] 



3.B. Adaptive Multichannel N on- Parametric Filters 

Adaptive Multichannel Non-Parametric Filters (AMNFs) [TT] approach the filtering problem 
from an estimation theoretic perspective. Specifically, these filters employ non-parametric 
kernel density estimators to determine the pixels in the filtered image as follows: 

y (r, c) - ^ x, ^ =ifc7 a* ((xo _ Xj)/fcj) ) 

n [a) 
hi = n K/3 £ || Xi-xj 1^ 
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where || . 1^ denotes the L\ (City-Block) norm. Two possible choices for the kernel function 
are the multivariate exponential K(x) = e - "*" 1 (AMNFE) and the multivariate gaussian 
K(x) = e" a5|ix|| 2 (AMNFG) functions. The scaling factor k in the kernel width calculation 
is set to the author recommended value of 0.33 [11]. The computational requirements of ([8]) 
can be reduced by speeding up the kernel computation. Both kernels involve the exponential 
(EXP) function which can be accurately approximated by polynomials. Note that in addition 
to AMNFs, the EXP function was used in the design of a number of other filters including the 
fuzzy vector median filter [23], fuzzy vector median-rational hybrid filter [19], kernel vector 
median filter [24], fast adaptive noise reduction filter [25], and self-adaptive noise reduction 
filter [26]. 

The argument of the EXP function in ([8]) depends on the k value, and the size and contents 
of a particular window. However, to obtain an accurate approximation, this argument needs 
to be constrained to a preferably small interval. Fortunately, for most practical purposes, we 
can set a cutoff point at T = 10.0 (e~ T = 4.539993e-05) and return for arguments outside 
the interval [0,T]. Figure H] shows a plot of the function in this interval. 
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Fig. 4. Function exp(-z) in the interval [0, 10] 



Table [2] shows the coefficients of the minimax polynomials of various degrees. Here, p and e 
represent the degree of the polynomial and the error of minimax approximation, respectively. 
It can be seen that the error values are relatively high, and as the approximation degree is 
increased, the accuracy doesn't improve significantly. This suggests that rational functions 
might be better suited for this approximation. Table |3] lists the coefficients of a minimax 
rational that approximates the EXP function with an error of e = 2.227050e-06. 
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Table 2. Minimax polynomials for the EXP function 



V 


£ 


ao 


Ol 


02 


0.3 


d4 


2 


1.785517e-01 


8.214528c-01 


-3.186948c-01 


2.544088c-02 






3 


8.259345e-02 


9.174126c-01 


-5.631179e-01 


1.015041c-01 


-5.519183c-03 




4 


3.337085e-02 


9.666313c-01 


-7.620584c-01 


2.145386c-01 


-2.509526c-02 


1.032877c-03 



Table 3. Minimax rational for the EXP function 



Term 


ao 


ai 


02 


03 


Ct4 


Numerator 


3.206619e-02 


-1.195191e-02 


1.756974c-03 


-1.199261e-04 


3.182685e-06 


Denominator 


3.206627e-02 


2.011147c-02 


5.853684e-03 


9.780143e-04 


1.251598e-04 



3. C. Entropy Vector Filters 

Entropy vector median filter (EVMF) introduced in [12] adaptively switches between the 
identity operation and a noise filtering mode to improve signal-detail preserving character- 
istics of standard filters such as the vector median filter, which performs fixed amount of 
smoothing in all pixel locations. Noise filtering is performed only in pixel locations which 
are identified as noisy by a switching operator. This is realized by comparing an adaptive 
threshold 0c expressed in the form of normalized entropy to a measure of normalized local 
contrast Pc as follows: 



y(7, c) 



Pi 



II*- 



argmm 

XiGW(r,c) 

x(r, c) 

x|| 2 




3 = 1 



\? 112 



'.Pi 



Pc>Pc 
otherwise 



(9) 



£?=i p 3 



where C = (n + 1 ) /2 and x denote the linear index of the center pixel (see Figure [T]) and 
the mean vector inside W(r, c), i.e. x = - Ya=i x *' respectively. 

Note that within the so-called generalized entropy vector filter (EVF) class [27], new filters 
can be designed by replacing the Euclidean distance function in §§§ with some other distance 
or similarity measure. 

The computational requirements of EVFs can be reduced by speeding up the entropy 
(ENT) function, whose argument falls into the interval [0,1]. Although, in theory, as the 
argument approaches 0, the function value approaches 0, in practice, this doesn't hold as 
the value of the logarithm function approaches negative infinity. Therefore, as in the case of 
the EXP function, we set a cutoff point at T = 0.05 and return for arguments less than T. 
Figure [5] shows a plot of the function in the interval [0.05, 1]. It can be seen that this function 
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is highly nonpolynomial [2], i.e. its derivatives are infinite at z = 0, and therefore using 
rational functions is more appropriate. Table H] lists the coefficients of a minimax rational 
that approximates the ENT function with an error of e — 7.342477e-07. 
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Fig. 5. Function zlog(z) in the interval [0.05, 1] 



Table 4. Minimax rational for the ENT function 



Term 


ao 


0.1 


0.2 


03 


04 


Numerator 


-1.519742c-04 


-6.835769c-02 


-8.856923c-01 


-5.369609c-01 


1.491165 


Denominator 


1.532270e-02 


3.987796c-01 


1.461793 


6.827004c-01 


-4.469776e-02 



4. Experimental Results 

In order to evaluate the performance and robustness of the presented approximations, a set 
of 100 high quality RGB images was collected from the Internet. The set included images of 
people, animals, plants, buildings, aerial maps, man-made objects, natural scenery, paintings, 
sketches, as well as scientific, biomedical, synthetic, and test images commonly used in the 
color image processing literature. 

The corruption in the test images was simulated using three noise models [28]: uncorrelated 
impulsive noise model, correlated impulsive noise model, and mixed noise model (Gaussian 




(10) 
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Noise + Correlated Impulsive Noise): 

Uncorrelated Impulsive Noise 
x = {x 1 ,x 2 ,x 3 } 

Ok with probability 1 — p^, , 
rfc with probability pk 
Correlated Impulsive Noise 

o with probability 1 — tp , 

{r 1; 02, 03} with probability (pi ■ tp, 
{01, r 2 , 03} with probability ip 2 • <p, 
{oi, o 2 , r%} with probability (p% ■ tp, 
{ri,r 2 , r 3 } with probability (1 - (ipx + <p 2 + ^3)) • V 9 

where o = {01,02,03} and x = {xi,x 2 ,Xs} represent the original and noisy color vectors, 
respectively, r = t 2 , r 3 } is a random vector that represents the impulsive noise, <p is 
the sample corruption probability, and ipi, p 2 , and ips are the corruption probabilities for 
the red, green, and blue channels, respectively. In the experiments, the channel corruption 
probabilities were set to 0.25. 

Filtering performance was evaluated by three effectiveness criteria j3]: 

M N 

1. Mean Absolute Error: MAE (X, Y) = 3^ £ II x ( r > c ) ~ y( r > c ) II 1 

r=l c=l 

where, X and Y denote respectively the M x N original and filtered images in the 
RGB color space. MAE measures the detail preservation capability of a filter. 

M N 

2. Mean Squared Error: MSE (X, Y) = 3^ £ II x ( r > c) - y(r, c) H2 

r=l c=l 

MSE measures the noise suppression capability of a filter. 



3. Normalized Color Difference: NCD(X,Y) 



M iV 

E E ll* I,ofc (r J c)-y £ «»(r,c)|| a 



E E ||x^(r,c)|| 2 

r— 1 c— 1 

where, x La6 (r, c) and y La6 (r, c) denote the CIEL*a*b* coordinates [1] of the pixel (r, c) 
in the original and filtered images, respectively. NCD measures the color preservation 
capability of a filter. 

The efficiency of a filter was measured by execution time in seconds (Programming Language: 
C, Compiler: gcc 3.4.4, CPU: Intel Pentium D 2.66Ghz). 

Table [5] shows the performance statistics for the three noise models. The test images 
were first corrupted using one of the noise models and then filtered using the exact and 
approximate versions of each filter. In the 'Mean' column, negative values and positive values 
for the MAE, MSE, and NCD indicate the percentage of filtering quality degradation and 
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improvement, respectively. For example, for 10% correlated impulsive noise, with respect to 
the MAE criterion, the approximate version of BVDF performs on the average 0.926% worse 
than the exact version, whereas with respect to the MSE criterion, the former performs 
0.171% better than the latter. On the other hand, for the execution time criterion, positive 
values indicate reduction in filtering time due to the use of the presented approximations. 
For example, the approximate version of BVDF is on the average 1371% (or 13.71 times) 
faster than the exact version. 



Table 5. Performance statistics at 10% noise level 





Uncorrelated Impulsive 


Correlated Impulsive 


Mixed 


Filter 


Measure 


Mean(%) 


Stdcv(%) 


Mean(%) 


Stdcv(%) 


Mean(%) 


Stdcv(%) 




MAE 


-1.000 


1.648 


-0.926 


1.558 


-0.022 


0.310 


BVDF 


MSE 


-0.571 


1.924 


0.171 


2.426 


-0.131 


1.103 


NCD 


-0.783 


1.283 


-0.759 


1.154 


0.007 


0.193 




Time 


1381.783 


16.576 


1371.314 


16.479 


1408.964 


11.964 




MAE 


-0.025 


0.137 


-0.016 


0.136 


-0.001 


0.137 


AMNFE 


MSE 


-0.070 


0.526 


-0.064 


0.754 


-0.006 


0.317 


NCD 


-0.020 


0.177 


-0.031 


0.195 


-0.004 


0.176 




Time 


142.700 


0.225 


143.496 


0.214 


140.814 


0.187 




MAE 


-0.141 


0.524 


-0.108 


0.489 


0.000 


0.133 


EVMF 


MSE 


-0.376 


1.732 


-0.203 


1.496 


0.013 


0.241 


NCD 


-0.147 


0.837 


-0.149 


0.625 


-0.004 


0.138 




Time 


236.582 


0.672 


236.105 


0.646 


215.000 


0.449 



It can be seen that in most cases the exact filters slightly outperform their respective 
approximate versions. This was expected since the approximate filters necessarily involve 
small amounts of computational error. Nevertheless, the difference between the approximate 
and exact versions for each filter is negligible for most practical purposes, which demonstrates 
the accuracy of the presented approximations. In addition, the low standard deviation values 
indicate the robustness of the approximations. 

The discrepancies in the speed up factors for the three filters can be attributed to the 
relative computational cost of the elementary functions involved. In other words, the speed 
up in BVDF is much greater than the other two filters because the ARCCOS function is 
computationally much more expensive than the EXP and ENT functions. 

Since the filters presented in Section [3] are primarily intended for the removal of impulsive 
noise, we conducted further experiments with the most commonly used impulsive noise 
model [2j[3], i.e. the correlated impulsive noise model [28]. Table |6] shows the performance 
statistics at 20%, 30%, and 40% noise levels. It can be seen that the performance of the 
approximate filters does not change significantly as the noise level is increased. 

Figure [6] compares the exact and approximate versions of each filter on the Lenna image. 
It can be seen that the presented approximations achieve substantial computational savings 
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Table 6. Performance statistics at higher noise levels 





20% 


30% 


40% 


Filter 


Measure 


Mean(%) 


Stdcv(%) 


Mean(%) 


Stdev(%) 


Mean(%) 


Stdcv(%) 




MAE 


-0.753 


1.311 


-0.626 


1.070 


-0.419 


0.861 


BVDF 


MSE 


0.193 


2.994 


0.027 


2.630 


-0.073 


1.855 


NCD 


-0.607 


1.016 


-0.530 


0.881 


-0.424 


0.802 




Time 


1342.074 


15.480 


1328.939 


14.784 


1301.070 


14.026 




MAE 


-0.032 


0.231 


0.030 


0.342 


-0.001 


0.315 


AMNFE 


MSE 


-0.125 


0.945 


0.021 


1.172 


0.006 


0.812 


NCD 


-0.045 


0.349 


0.034 


0.478 


0.027 


0.395 




Time 


143.460 


0.241 


145.513 


0.226 


145.362 


0.249 




MAE 


-0.030 


0.676 


0.013 


0.534 


0.067 


0.476 


EVMF 


MSE 


-0.056 


1.539 


0.046 


1.105 


0.065 


0.867 


NCD 


-0.004 


0.829 


0.026 


0.633 


0.088 


0.547 




Time 


228.346 


0.551 


222.356 


0.431 


216.354 


0.423 



without introducing any perceivable artifacts on the filtering results. In addition, the MAE 
and MSE values indicate that the filtering effectiveness of the exact and approximate filters 
are virtually the same. 

5. Conclusions 

In this article, we proposed a novel approach to speed up popular vector filters using minimax 
approximations. Advantages of this approach include ease of implementation, extremely good 
accuracy, and high computational speed. The presented approach can be adapted to other 
noise removal filters that involve computationally expensive mathematical functions. Finally, 
the given approximations have applications that go beyond color image filtering including 
computer graphics and computational geometry 

Acknowledgments 

This publication was made possible by a grant from The Louisiana Board of Regents 
(LEQSF2008-11-RD-A-12). 

References 

1. K.N. Plataniotis and A.N. Venetsanopoulos, Color Image Processing and Applications 
(Springer, 2000). 

2. R. Lukac and K.N. Plataniotis, "A Taxonomy of Color Image Filtering and Enhancement 
Solutions," in Advances in Imaging & Electron Physics (Vol. 140), P.W. Hawkes, ed. 
(Academic Press, 2006), pp. 187-264. 



Journal of the Optical Society of America A, 26(6): 1518-1524, 2009 




(a) Original (512 x 512) (b) 10% noisy (MAE 6.373, (c) BVDF exact (MAE 

MSE 987.418, Time 0) 3.936, MSE 43.571, Time 

10.360) 




(d) BVDF approx (MAE (e) AMNFE exact (MAE (f) AMNFE approx (MAE 
3.936, MSE 43.558, Time 3.471, MSE 30.303, Time 3.471, MSE 30.303, Time 
0.688) 0.594) 0.422) 




(g) EVMF exact (MAE (h) EVMF approx (MAE 
1.139, MSE 25.898, Time 1.138, MSE 25.819, Time 
0.594) 0.250) 



Fig. 6. Comparison of the exact and approximate filters on the Lenna image 



3. M.E. Celebi, H.A. Kingravi, and Y.A. Aslandogan, "Nonlinear Vector Filtering for Im- 
pulsive Noise Removal from Color Images," Journal of Electronic Imaging 16(3), 033008 



Journal of the Optical Society of America A, 26(6): 1518-1524, 2009 

(2007). 

4. R. Lukac, B. Smolka, K. Martin, K.N. Plataniotis, and Venetsanopoulos A.N., "Vector 
Filtering for Color Imaging," IEEE Signal Processing Magazine 22(1), 74-86 (2005). 

5. M.E. Celebi and Y.A. Aslandogan, "Robust Switching Vector Median Filter for Impulsive 
Noise Removal," Journal of Electronic Imaging 17(4), 043006 (2008). 

6. M. Barni, "A Fast Algorithm for 1-Norm Vector Median Filtering," IEEE Trans, on 
Image Processing 6(10), 1452-1455 (1997). 

7. M. Barni, F. Buti, F. Bartolini, and V. Cappellini, "A Quasi-Euclidean Norm to Speed 
up Vector Median Filtering," IEEE Trans, on Image Processing 9(10), 1704-1709 (2000). 

8. J. Astola, P. Haavisto, and Y. Neuvo, "Vector Median Filters," Proc. of the IEEE 78(4), 
678-689 (1990). 

9. M.E. Celebi, H.A. Kingravi, B. Uddin, and Y.A. Aslandogan, "Fast Switching Filter for 
Impulsive Noise Removal from Color Images," Journal of Imaging Science and Technol- 
ogy 51(2): 155-165 (2007). 

10. P.E. Trahanias and A.N. Venetsanopoulos, "Vector Directional Filters: A New Class of 
Multichannel Image Processing Filters," IEEE Trans, on Image Processing 2(4), 528-534 
(1993). 

11. K.N. Plataniotis, D. Androutsos, S. Vinayagamoorthy, and A.N. Venetsanopoulos, 
"Color Image Processing Using Adaptive Multichannel Filters," IEEE Trans, on Image 
Processing 6(7), 933-949 (1997). 

12. R. Lukac, B. Smolka, K.N. Plataniotis, and A.N. Venetsanopoulos, "Entropy Vector 
Median Filter," in Proc. of the IbPRIA Conf., Lecture Notes in Computer Science 2652, 
pp. 1117-1125 (2003). 

13. E.W. Cheney, Introduction to Approximation Theory (AMS, 2000). 

14. J.-M. Muller, Elementary Functions: Algorithms and Implementation (Birkhauser, 2006). 

15. W. Fraser, "A Survey of Methods of Computing Minimax and Near-Minimax Polynomial 
Approximations for Functions of a Single Independent Variable," Journal of the ACM 
12(3), 295-314 (1965). 

16. N. Nikolaidis and I. Pitas, "Nonlinear Processing and Analysis of Angular Signals," IEEE 
Trans, on Signal Processing, 46(12), 3181-3194 (1998). 

17. P.E. Trahanias, D. Karakos, and A.N. Venetsanopoulos, "Directional Processing of Color 
Images: Theory and Experimental Results," IEEE Trans, on Image Processing 5(6), 868- 
880 (1996). 

18. D.G. Karakos and P.E. Trahanias, "Generalized Multichannel Image Filtering Struc- 
tures," IEEE Trans. Image Processing 6(7), 1038-1045 (1997). 

19. L. Khriji and M. Gabbouj, "Adaptive Fuzzy Order Statistics- Rational Hybrid Filters for 
Color Image Processing," Fuzzy Sets and Systems 128(1), 35-46 (2002). 



Journal of the Optical Society of America A, 26(6): 1518-1524, 2009 

20. R. Lukac, "Adaptive Color Image Filtering Based on Center- Weighted Vector Directional 
Filters," Multidimensional Systems and Signal Processing 15(2), 169-196 (2004). 

21. K.N. Plataniotis, D. Androutsos, and A.N. Venetsanopoulos, "Color Image Processing 
Using Adaptive Vector Directional Filters," IEEE Trans, on Circuits and Systems-II 
45(10), 1414-1419 (1998). 

22. R. Lukac, B. Smolka, K.N. Plataniotis, and A.N. Venetsanopoulos, "Vector Sigma Filters 
for Noise Detection and Removal in Color Images," Journal of Visual Communication 
and Image Representation 17(1), 1-26 (2006). 

23. K.N. Plataniotis, D. Androutsos, and A.N. Venetsanopoulos, "Adaptive Fuzzy Systems 
for Multichannel Signal Processing," Proceedings of the IEEE 87(9), 1601-1622 (1999). 

24. B. Smolka, K.N. Plataniotis, R. Lukac, and A.N. Venetsanopoulos, "Kernel Density 
Estimation Based Impulsive Noise Reduction Filter," in Proc. of the IEEE Int. Conf. on 
Image Processing (ICIP'03) 2, pp. 137-140 (2003). 

25. B. Smolka, R. Lukac, A. Chydzinski, K.N. Plataniotis, and K. Wojciechowski, "Fast 
Adaptive Similarity Based Impulsive Noise Reduction Filter," Real-Time Imaging 9(4), 
261-276 (2003). 

26. B. Smolka, K.N. Plataniotis, A. Chydzinski, M. Szczepanski, A.N. Venetsanopoulos , 
K. Wojciechowski, "Self-Adaptive Algorithm of Impulsive Noise Reduction in Color Im- 
ages," Pattern Recognition 35(8), 1771-1784 (2002). 

27. R. Lukac, B. Smolka, K.N. Plataniotis, and A.N. Venetsanopoulos, "Generalized Entropy 
Vector Filters," in Proc. of the 4-th EURASIP EC-VIP-MC, Video, Image Processing and 
Multimedia Communications Conf., pp. 239-244 (2003). 

28. T. Viero, K. Oistamo, and Y. Neuvo, "Three- Dimensional Median-Related Filters for 
Color Image Sequence Filtering," IEEE Trans, on Circuits and Systems for Video Tech- 
nology 4(2), 129-142 (1994). 



